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(54) System and method for mapping chromatic dispersion in optical fibers 



(57) A system and method in which the chromatic 
dispersion characteristic, as a junction of distance along 
a section of transmissive fiber, is measured. First and 
second optical signals are launched into a fiber under 
test to thereby generate, by a tour-wave mixing process 
in the fiber, a probe signal. Because of a wave-vector 
phase mismatch, the probe signal power oscillates with 
a spatial frequency that can be measured as a function 



of distance in the fiber. The intensity oscillations are 
measurable as, for example, temporal variations in the 
Rayleigh backscattered light detected at the input end 
of the fiber. According to the present invention, the dis- 
persion parameter at one or both of the first and second 
optical signal wavelengths, as a function of length along 
the fiber, is derived directly from these intensity oscilla- 
tion measurements. From this information, it is possible 
to further derive maps at other wavelengths of interest. 
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Description 

Field Of The Invention 

5 The present invention relates generally to optical fiber communications and, more particularly, to systems and 

methods for measuring and managing chromatic dispersion in optical fibers. 

Background Of The Invenlion 

10 Fiber chromatic dispersion has played an important role in the design of optical fiber systems for more than a 

decade. Until the advent of the erbium-doped fiber amplifier, the systems were more or less linear. Hence, it was only 
the integrated dispersion over a fiber span that influenced system performance. As the need to satisfy the demand for 
transmission capacity over greater distances has led to more sophisticated processing of optical signals, dispersion 
management - in which a dispersion "map" is chosen to minimize and/or harness the effects of fiber nonltnearities, 

15 has becoming an increasingly important tool. 

In dispersion -shifted fiber (DSF), dispersion is known to vary as a function of location in the fiber. In a paper by K. 
Inoue entitled "Four-Wave Mixing in an Optical Fiber in the Zero-Dispersion Wavelength Region", J. Lightwave Technol., 
Vol. 10, pp. 1553-1561 (1992), for example, it was reported that when a 10-km section of DSF was cut intofour 2.5-km 
segments, the average dispersion zero wavelength for the segments for the segments varied by at least 1 nm - a 

20 significant deviation for some applications. Accordingly a reliable map of chromatic dispersion can not be obtained 
merely by measuring the average dispersion in the fiber span. 

It has therefore been proposed to map the distribution of chromatic dispersion along a fiber span using a Rayleigh 
backscattering technique. This destructive technique, which relies on the dependence of the dispersion zero on the 
fiber core size, is described in an article by M. Ohashi and M. Tateda entitled "Novel Technique for Measuring the 

25 Distributed Zero- Diapers ion Wavelength of Optical Fibers", Electron Letters., 29, 426-428 (1993). If the doping of the 
fiber preform does not change over its length, then changes in the dispersion zero can be inferred from changes in the 
core size. Changes in core size are estimated using optical time domain reflectometry (OTDR) to determine the capture 
ratio for Rayleigh backscattered light. By summing OTDR measurements taken in opposite directions, the effects of 
fiber loss are removed and changes in the capture ratio are observed and used to determine variations in the fiber 

30 dispersion zero. 

Recently, a non-destructive dispersion measurement method was described that determined the local dispersion 
zero from modulation-instability-induced gain at wavelengths longer than the dispersion zero. A strong pump pulse of 
wavelength and a weak signal pulse of wavelength are injected simultaneously into a test fiber with the difference 
between >.p and being about 5 to 10 nm. The backscattered signal light is observed through OTDR. When the pump 

35 wavelength is near the dispersion zero, but in the anomalous dispersion region, the modulation instability will provide 
gain for the probe pulse - gain that can be observed in the OTDR trace. Thus, reduction in the slope of the OTDR at 
a particular distance into the fiber indicates that the pump is experiencing anomalous dispersion at that point in the 
fiber. To map fiber dispersion, the pump and probe wavelength are swept, maintaining a constant separation, Xp - X^. 
and the resulting OTDR traces are recorded. The dispersion zero of a particular point in the fiber is at the short- 

40 wavelength side ot those pump wavelengths for which modulation-instability gain is observed. This technique has 
demonstrated a wavelength resolution of 0.2 nm and a spatial resolution of about a kilometer 

Yet another technique that has been proposed uses partially-degenerate four-photon mixing to determine the dis- 
persion zero. Essentially, the mixing generates an idler wave from pump and signal waves of angular frequencies oOp 
and cOg propagating in the fiber. The power ot the idler wave with frequency o)i=2o)p - cOg, will be maximized when the 

^5 process is phase matched, that is, when Ap= 2|3(a)p) - p(cos) * P(^) = 0, where p(o)) is the propagation constant. To first 
approximation, phase matching occurs when cOp is set to the dispersion zero of the fiber Thus, by tuningtOp, and looking 
for a maximum in idler power, it is possible to measure the zero-dispersion wavelength. Distance resolution is obtained 
by using signal and pump pulses with widely-separated wavelengths. Specifically, the differing group velocities of the 
pump and signal pulses cause the pump to overtake the signal pulse (assuming the pump wavelength is near the 

so dispersion zero of the fiber). With sufficient group-velocity dispersion and short enough pulses, the region of overlap 
of the pulses occurs over some useful distance. The timing of the pulses at a>p and can then be adjusted so that 
this overlap occurs at some desired point within the fiber 

Each of the above-described techniques penmits measurement of the wavelength of zero dispersion. Disadvan- 
tageously: such techniques require extensive data-gathering over a considerable wavelength range, so that the meas- 

55 urements for just one dispersion map take a long time. Moreover, access is required to both ends of the fiber under 
investigation. 
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Method and apparatus according to the invention are as set out in the independent claims. Preferred forms of the 
method are set out in the dependent claims, preferred forms of the apparatus being adapted to carry out the preferred 
5 forms of the method. 

The aforementioned deficiencies are addressed, and an advance is made in the art, by a system and method in 
which the chromatic dispersion parameter, as a function of distance along a section of transmissive fiber, is measured. 
In accordance with an illustrative embodiment of the inventbn first and second optical pulses are repetitively launched 
into a fiber under test to thereby generate, by a four-wave mixing process in the fiber, a probe signal. Because of the 

10 wave-vector phase mismatch, the probe signal power oscillates with a spatial frequency that can be measured as a 
function of distance in the fiber. These intensity oscillations are measurable as, for example, as temporal variations in 
the Rayleigh backscattered light detected at the input end of the fiber According to the present invention, the dispersion 
parameter at one or both of the first and second optical signal wavelengths, as a function of length along the fiber, is 
derived directly from these intensity oscillation measurements. From this information, it is possible to further derive 

IS maps at other wavelengths of interest. 



Brief Description of the Drawings 

Various advantages of the present invention will become apparent to one skilled in the art upon by reference to 
20 the specification which follows and to drawings, in which: 

FIG. 1 is a graph showing the spectrum of the four-wave mixing process employed by the present invention; 

FIG. 2 is a graph depicting the most probable band for transmission, and most suitable wavelengths for measure- 
rs ment of dispersion parameter maps, in relation to the gain profile of an erbium fiber amplifier; 

FIG. 3 is a block diagram showing a chromatic dispersion characteristic measuring system according to an illus- 
trative embodiment of the present invention. 

20 FIG. 4 is a graph depicting several samples of the signal retumed from a "dispersion tapered" span of fiber; and 

FIG. 5 is a graph comparing dispersion maps obtained utilizing the method of the present invention with those 
obtained by a conventional, destructive technique. 



25 Detailed Description of the Invention 



According to the present invention, strong, sub -microsecond pulses at fixed wavelengths X-^ and respectively, 
are simultaneously launched into a length of fiber under test, so that they may generate respective FWM product fields 
at the Stokes wavelength and the anti-Stokes wavelength Xp^, sequentially in each part of the fiber. By measuring 
40 the frequency of the intensity oscillations in the probe signals generated in this manner, a detailed dispersion map, D 
(Xj,z) can be obtained directly, where D is the dispersion parameter, X-^ is the wavelength of one of the two aforemen- 
tioned fixed-wavelength sources used and z is the distance to a particular point along the fiber Before proceeding to 
a description to an illustrative system which may be utilized to acquire a detailed dispersion map in accordance with 
the present invention, the fundamental principles involved will first be discussed briefly. 

45 

Chromatic Dispersion 

The dispersion relation for an optical fiber is given by: 



so 



55 
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(00 -CDo) + l/2 



(a)-a)(>)^+. 



(1) 



where ko and all derivatives are evaluated at the (arbitrary) frequency o]^. The first derivative of the dispersion relation 
with respect to to is the inverse group velocity, and its second derivative is the corresponding dispersion. The dispersion 
parameter. D, is the wavelength derivative of the inverse group velocity: 
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5 Four-Wave Mixing Processes 



Essentially, the present invention is based upon the realization by the inventors herein that the dispersion parameter 
D(Xi) Is directly proportional to the wave-vector nnisnnatch for one of two four-wave mixing processes depicted in FIG. 
Una first of these processes, two photons at angular frequency cOj combine with one at angular frequency to form 
10 a Stokes photon at cOg. Symmetrically, in a second mixing process two photons at angular frequency cog combine with 
one at angular frequency a>i to form an anti-Stokes photon at co^. Thus, for energy conservation, all spectral terms in 
FIG. 1 are separated uniformly by an amount 5co. 

For a detailed discussion of the well known physics of four-wave mixing, reference may be had to a paper by the 
inventors herein entitled "Pseudo-phase-matched four-wave mixing in soliton WDM transmission", Opt Lett 21, 396 
15 (1996). 



Wave-Vector Mismatch 



It can then be shown that the wave-vector mismatch. 5/c, for the aforementioned first and second mixing process 
20 is in direct proportion to dispersion parameters D{X^) and Di)^), respectively In the case of the first mixing process, 
for example, the wave-vector mismatch can be related to the dispersion parameter D(^-,) as foltows: 



25 



5©' ='2ncD{X,)\ 



(3) 



The final form of Eq. (3) is accurate because in accordance with the inventive method, the condition 5oV(o = -SX/X « 
1 is always satisfied. It should also be noted that the net term in Eq. (3) contributed by third order dispersion d^k/aco^ 
oc dD/dX is identically zero. The effects of still higher orders of dispersion, while not always identically zero, nevertheless 
tend to be negligible because of the small sizes of the derivatives and of the interval 6X. 



Calculation of the Signal Power 

As Indicated above, the inventive technique of the present invention involves launching two optical signals to 
thereby generate respective FWM product fields at the Stokes wavelength X^ and the anti-Stokes wavelength Xp^ se- 
quentially in each part of the fiber. Because of the aforementioned phase mismatch (bk), the corresponding power of 
the Stokes signal Ps(z) oscillates with the spatial frequency 

F3=1/A3=g=cD(X,)(f)^ (4) 



Accordingly, by measuring the frequency Fs(z) of the intensity oscillations in the power of the Stokes signal, It is possible 
45 to Obtain a dispersion map 0(X^,z) having a spatial resolution As. Of course, it will be readily appreciated by those 
skilled in the art that it is equally possible to obtain a dispersion map D{X2,z) by measuring the frequency Fa(z) of 
intensity oscillations in the power of the antl-Stokes signal. As such, an accurate dispersion map for one or both of the 
signals launched into the fiber under test can thus be readily obtained by applying the teachings of the present invention. 
Moreover, from knowledge of the third order dispersion constant - which can be quickly determined by comparing 
so measurements of dispersion maps D(X^,z) and D(X2;Z) -- either of the measured dispersion maps can be easily con- 
verted to that for any other desired wavelengths. 

A variety of techniques may be utilized to measure the aforementioned intensity oscillations. By way of Illustrative 
example, they may be observed in Rayleigh backscattering at the input end of the fiber under test. There, the signal 
will fluctuate at a temporal frequency 

55 

^sigO^ = fnW 
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where n is the effective index of refraction of the fiber, and where t is the roundtrip time from the fiber input to point 2 
and relum, i.e.. 

t-'-^ (6, 

For a typical index of 1 .46, St is 9.73 jis for each kilometer of fiber Combining Eqs. (4) and (5). the dispersion map is 
obtained for the first input signal at wavelength 

0(.,.z)=y(J-^rW;=^"z) (7, 

As will be discussed in detail later, the frequencies dictated by Eq. 7 are typically in a range of some tens to a few 
IS hundreds of kilohertz. 

From the known loss and scattering properties of the fiber, and an understanding of the physics of the aforemen- 
tioned four-wave mixing process, the strength of the Rayleigh backscattered signal can be estimated. For the case 
where the input pulses at and are co-polarized, and where there is no significant initial signal at and X^^, the 
signal power at the Stokes wavelength X3 is 



10 



20 



o^2 

(8) 



25 As will be immediately apparent to those skilled in that art, a similar expression yields the signal power P>v(z), at the 
anti-Stokes wavelength X^^. In Eq. (8), P^^ and are the pulse powers at the fiber input, is the effective area of 
the fiber core, is the non-linear index coefficient (n^ is 2.7 X 10"^ |j.m2/W for silica glass fibers), R \s the Rayleigh 
back-scattering coefficient, bz is the fiber length occupied by the pulses at any given time, and a is the fiber's loss 
coefficient. The factor of four in the exponential loss term stems from the facts that (1) the product P'i^P2 declines as 

20 exp(-3otz). and (2) the Rayleigh backscattering at a given point z along the fiber suffers an additional loss factor of exp 
i-az) in retuming to the fiber input. From Eq. (8). it can be seen that for pulse input powers {P^^ and P2°) on the order 
of 1 W, there should be adequate signal strength for measurements of fiber spans up to several tens of kilometers long. 

With particular reference now to FIG. 2. it will be appreciated that the band of wavelengths currently of greatest 
interest for long distance transmission corresponds to the broad, relatively flat gain peak of erbium fiber amplifiers lying 

35 roughly between 1552 and 1560 nm. For soliton transmission using "dispersion tapered" fiber spans, the wavelength 
of zero dispersion, Xq, will correspondingly lie in the range between --1542 to 1553 nm, and for "non-retum-to-zero 
systems" (NRZ), X^ will lie at even longer wavelengths. For several reasons, it is preferred, though not required, that 
the input or measuring wavelengths X^ and X2 be to the short wavelength side of the shortest anticipated value for Xq. 
This relationship will, for example, ensure that there is no "modulational instability", regardless of the respective power 

40 levels. At power levels approaching 1 W, this would certainly not be the case if X•^ and were greater than X^. Moreover, 
since the signals do not depend on the sign of the dispersion parameter D, the dispersion map tends to be ambiguous 
unless it is always well to one side of D=0. Finally, the spatial resolution tends to become poor as D(A^) approaches 
zero. The inventors herein have observed that excellent results are obtained, for example, when X^ and are chosen 
to lie in the vicinity of the 1530 nm gain peak of an erbium doped amplifier, as shown in FIG. 2. Of course, if the 

45 wavelength band of interest is not to the short wavelength side of the shortest anticipated value forXo, as it would not 
be in the popular case where the band of >^ is centered about, say 1515 nm. then the input power levels of the measuring 
wavelength signals should be reduced below the threshold for modulational instability, that is, below about 200 mw. 
Alternatively, an amplifier with a different characteristic gain curve might be employed. 

It should be noted that there is a non-linear contribution to the wave-vector mismatch dk, which for the first (Stokes) 

50 mixing process is given by the relation 

5/r„,=Y(2P, .P2) (9) 

with a similar expression (just reverse the subscripts 1 and 2) for the second (anti-Stokes) mixing process. In Eq. (9), 
y is given by the relation 
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2r, "2 



off 



(10) 



5 For the first mixing process. S/c^i is zero if Pg = 2?^. Identically, of course. S/c^i is zero if P-, = 2P2 for the second mixing 
process. Even if the appropriate one of those conditions is not precisely met, however, note that bk^^ will tend to be 
only a small fraction of the linear bkior all dispersion parameters D but those very close to zero, since y is numerically 
very small about 

2.5 W'^km-i where A^fj = 50 iim^ and his approximately 1530 nm, for example. 

TO With reference now to FIG. 3, there is shown a dispersion map measuring system 10 constructed in accordance 

with an illustrative embodiment of the present invention. The fixed input signals and 7^ are generated by first and 
second cw sources 12 and 14, respectively In the illustrative embodiment, the sources 12 and 14 are configured as 
diode lasers tuned to provide pulses at 1532.6 nm and 1536.15 nm, respectively (for a wavelength separation of 3.55 
nm). It will, of course, be readily appreciated by those skilled in the art that the wavelengths provided are for exemplary 

ts purposes only, and, moreover, they may be generated by any suitable optical signal generating device, including an 
optical pulse generation. Because orthogonal polarizations have the potential to reduce the signal intensity by up to a 
factor of ten or more, care should be taken to see that the cw sources 1 2 and 1 4 are co-polarized to provide an adequate 
signal intensity. 

The cw signals output by optical sources 12 and 14 are combined and input to a phase modulator 16. While not 

20 required, the phase modulator advantageously serves to remove the annoying high-frequency 1 MHz) modulation 
from the input signals that results from the finite temporal width of the exciting pulses. The output of the phase modulator 
is amplified by amplifier 18, and an accousto-optic modulator 20 coupled to a pulse generator 22 shapes the output of 
the cw sources into pulses at and X^. The output of the accousto-optic modulator 20 is supplied to an optical amplifier 
24, an erbium doped fiber amplifier (EDFA) in the illustrative example, which amplifier provides sufficient gain to bring 

25 the pulse signals up to peak power levels near 1 W. Even though the pulse powers at the output of amplifier 24 are on 
the order of 1 W, the low duty cycle and energy storage performance of the amplifier allow a pump power for the 
amplifier that is on the order of only a few tens of mw. 

A three-port optical circulator 26 launches the pulses into the fiber under test, indicated generally at F, and directs 
the Rayleigh backscattered signal to a tunable receiver arrangement 28. In the illustrative apparatus depicted in FIG. 

30 3, the receiver arrangement 28 comprises an EDFA pre-amplifier 30 which provides sufficient gain to allow subsequent 
processing of the backscattered signals. Tuning to the desired backscattered signal, which will at either the Stokes or 
anti-Stokes wavelength, is achieved by a tunable narrow band optical filter 32 which passes the signal of interest to a 
photodetector 34 for detection. Such an optical filter may be constructed, for example, by cascading a plezo-tuned 
Fabry-Perot etaion with a bandwidth of 20 Ghz and a free spectral range of 16 nm, with a tuneable, 1 nm bandwidth 

35 interference filter (neither of which are shown). The output of photodetector 34 is filtered by a bandpass filter 36, 
amplified by an amplifier 37, and supplied to a sampling osclliiscope 38. 

FIG. 4 shows samples of the signal returned from a 34.4 km long, "dispersion tapered" span of fiber under test 
using the apparatus of FIG. 3. The samples were of the signal voltage waveform at = 1 539.70 nm. and as averaged 
over a few hundred repetitions with a required total averaging time of less than 0.2 seconds. The high-intensity input 

40 end is at 0 km. An excellent signal-to-noise ratio can be observed, even at the far end of the span. This is due, in part, 
to the narrow effective bandwidth used to compute the noise referred to the input of the EDFA pre-amplifier 30 (FIG. 
3), which was about 70 Mhz taken from the geometric mean of the bandwidths, 20 Ghz and 250 kHz, respectively of 
the two filters in the receiver chain. Similar results were obtained for the anti-Stokes signal waveform at = 1 529.05 nm. 
With reference now to FIG. 5, there are shown the nearly identical dispersion maps D(1557, z) of the same dis- 

45 perston tapered span discussed above, as inferred from (1) the Stokes signal, (2) the ant-Stokes signals, and (3) as 
translated to the common wavelength of 1 557 nm. The dashed curve corresponds to the ideal exponential taper. Similar 
results were also obtained using a conventional destructive technique during which the fiber was cut into five sections. 
The results obtained in accordance with the present invention compare favorably to the results of the destructive ex- 
amination, the latter being shown by the solid horizontal lines in FIG. 5. For the translation, third order dispersion values 

so were used, as determined from a comparison of the Stokes and anti-Stokes signals themselves. These third order 
dispersion values were not strictly constant across the span. 

For the determination of the dispersion maps shown In FIG. 5, the pump intensities at and X2 were adjusted to 
be in a two-to-one ratio, respectively, for the measurement at X^, and in a one-to-two ratio for the measurement at Xf^. 
According to Eq. (9), supra, no nonlinear contribution to the wave-vector mismatch should be discemible. This prediction 

55 was confirmed by noting a lack of any measurable change in the observed frequencies as the net pump power is 
reduced, so long as the two-to-one power ratios are maintained. 

In view of the foregoing, it should be readily apparent to those skilled in the art that the inventive technique described 
herein makes it possible to obtain rapid and facile measurements of the dispersion maps of long fiber spans, with an 
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40 



accuracy of +/-0.03 ps/nm-km or better throughout the span, and with a spatial resolution of a fraction of a kilometer 
Although one specific embodiment of the present invention has been described, various changes and modifications 
may be made by those skilled in the art without departing from the scope of the present invention, as defined by the 
appended claims. 

Claims 

1, A method of mapping chromatic dispersion in a span of optical fiber, comprising the steps of: 



launching Into a first end of said fiber, simultaneously, a first optical signal at a first wavelength and a second 
optical signal at a second wavelength to thereby generate, via a four wave mixing process in which two photons 
at said first wavelength combine with one at said second wavelength, a probe signal; 

measuring the frequency of the intensity oscillations, as a function of distance in the fiber, in a sample of the 
?5 probe signal; and 

deriving, for at least one of the first and second wavelengths, a chromatic dispersion parameter as a function 
of distance along the fiber from the measurement obtained during the measuring step. 

2. The method of claim 1 , wherein said probe signal is a Stokes signal, the second wavelength is shorter than the 
20 first wavelength and the ratio of launched power between the first and second wavelengths is approximately 1 :2. 

3. The method of claim 1 , wherein said probe signal is an anti-Stokes signal, the first wavelength is shorter than the 
first wavelength and the ratio of launched power between the first and second wavelengths is approximately 2:1 . 

25 4. The method of any of the preceding claims, wherein a chromatic dispersion parameter is derived for each of said 
first and second wavelengths and wherein said method further includes a step of calculating a third order dispersion 
parameter from differences between the two derived chromatic dispersion parameters. 

5. The method of any of the preceding claims, wherein said measuring step is performed by observing a Rayleigh- 
20 backscattered sample of the probe signal at the first end of the fiber. 

6. The method of any of claims I to 4, wherein said measuring step is performed by observing a plurality of side- 
scattered samples of the probe signal. 

2S 7. The method of any of the preceding claims, in which the first and second wavelengths are shorter than the shortest 
anticipated wavelength of zero dispersion in the span under test. 

8. The method of any of the preceding claims, wherein said launching step comprises repetitively launching pulses 
at said first and second wavelengths. 



9. An apparatus for measuring the chromatic dispersion characteristic in a span of optical fiber, comprising: 



an optical signal generating arrangement for launching into a first end of said fiber, simultaneously, a first 
optical signal at a first wavelength and a second optical signal at a second wavelength to thereby generate, 
45 via a four wave mixing process in which two photons at said first wavelength combine with one at said second 

wavelength, a probe signal; and 

a detecting arrangement for measuring the frequency of the intensity oscillations, as a function of distance in 
the fiber, in a sample of the probe signal, whereby for at least one of the first and second wavelengths, a 
chromatic dispersion parameter can be derived as a function of distance along the fiber. 

so 
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FIG. 1 
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(54) System and method for mapping chromatic dispersion in optical fibers 



(57) A system and method In which the chromatic 
dispersion characteristic, as a function of distance along 
a section of transmissive fiber, is measured. First and 
second optica) signals are launched into a fiber under 
test to thereby generate, by a four-wave mixing process 
in the fiber, a probe signal. Because of a wave-vector 
phase mismatch, the probe signal power oscillates with 
a spatial frequency that can be measured as a function 



of distance in the fiber. The intensity oscillations are 
measurable as, for example, temporal variations in the 
Rayleigh backscattered light detected at the input end 
of the fiber. According to the present invention, the dis- 
persion parameter at one or both of the first and second 
optical signal wavelengths, as a function of length along 
the fiber, is derived directly from these intensity oscilla- 
tion measurements. From this information, it is possible 
to further derive maps at other wavelengths of interest. 



FIG. 3 
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